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Measuring the proper motions and geometric distances of galaxies within the Local Group is 
very important for our understanding of the history, present state and future of the Local Group. 
Currently, proper motion measurements using optical methods are limited only to the closest com- 
panions of the Milky Way. However, Very Long Baseline Interferometry (VLBI) provides the best 
angular resolution in astronomy and phase-referencing techniques yield astrometric accuracies of 
« 10 micro-arcseconds. This makes a measurement of proper motions and angular rotation rates 
of galaxies out to a distance of ~ 1 Mpc feasible. 

This article presents results of VLBI observations of regions of H2O maser activity in the Local 
Group galaxies M33 and IC 10. Two masing regions in M33 are on opposite sides of the galaxy. 
This allows a comparison of the angular rotation rate (as measured by the VLBI observations) 
with the known inclination and rotation speed of the HI gas disk. This gives a geometric distance 
of 730 ± 100 ±135 kpc. The first error indicates the statistical error from the proper motion mea- 
surements while the second error is the systematic error from the rotation model. This distance 
is consistent, within the errors, with the most recent Cepheid distance to M33. Since all position 
measurements were made relative to an extragalactic background source, the proper motion of 
M33 has also been measured. This provides a three dimensional velocity vector of M33, showing 
that this galaxy is moving with a velocity of 190 ± 59 km s^' relative to the Milky Way. For 
IC 10, we obtain a motion of 215 ±42 km s^' relative to the Milky Way. These measurements 
promise a new handle on dynamical models for the Local Group and the mass and dark matter 
halo of Andromeda and the Milky Way. 
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1. Introduction 

The nature of spiral nebulae like M33 was the topic of a famous debate in the 1920's. While 
some astronomers favored a close distance and Galactic origin, others were convinced of the ex- 
tragalactic nature. In 1923, van Maanen claimed to have measured a large proper motion and 
angular rotation of M33 from photographic plates separated by 12 years [|l]]. These measure- 
ments yielded rotational motions of ^ 10-30 mas yr^^, clearly indicating a close distance to M33. 
However, Hubble discovered Cepheids in M33, showing a large distance to M33 and confirming 
that M33 is indeed an extragalactic object The expected proper motions from the rotation of 
M33 are only f» 30 /las yr^\ 3 orders of magnitude smaller than the motions claimed by van 
Maanen. The source of error in van Maanens motion determinations was never identified. After 
more than 80 years, the idea behind the experiment to measure the rotation and proper motions 
of galaxies remains interesting for our understanding of the dynamics and geometry of the Local 
Group. Hence, they are an important science goal of future astrometric missions (e.g. Gaia [||]). 

The problem when trying to derive the gravitational potential of the Local Group is that usu- 
ally only radial velocities are known and hence statistical approaches have to be used. Kulessa 
and Lynden-Bell introduced a maximum Ukelihood method which requires only the line-of-sight 
velocities, but it is also based on some assumptions (eccentricities, equipartition) [Q]. 

Clearly, the most reliable way of deriving masses is using orbits, which requires the knowledge 
of three-dimensional velocity vectors obtained from measurements of proper motions. However, 
measuring proper motions of members of the Local Group to determine its mass is difficult. For the 
LMC aproper motion of 1.2ib0.3 mas yr^' was obtained from comparing photographic plates over 
a time-span of 14 years The Sculptor dwarf spheroidal galaxy moves with 0.56 ±0.25 mas yr^' 
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obtained from plates spanning 50 years in time (Fig |l|, [^). It was shown that the inclusion of these 
marginal proper motions can already significantly improve the estimate for the mass of the Milky 
Way, since it reduces the strong ambiguity caused by Leo I, which can be treated as either bound 
or unbound to the Milky Way [^]. 

In recent years, the proper motions of a number of Galactic satellites have been measured 



using the HST 081 R, H^, 11, 12, 13, 14, hj, 16 1. These galaxies are all closer than 150 kpc and 



show motions between 0.2 and a few milliarcseconds (mas) per year. More distant galaxies, such 
as galaxies in the Andromeda subgroup at distances of ^ 800 kpc, have smaller angular motions, 
which are currently not measurable with optical telescopes. 

2. Proper Motions with VLBI 

On the other hand, the expected proper motions for galaxies within the Local Group, ranging 
from 0.02 - 1 mas yr^^, can be seen with Very Long Baseline Interferometry (VLBI) using the 
phase-referencing technique. A good reference point is the motion of Sgr A* across the sky reflect- 
ing the Sun's rotation around the Galactic Center. This motion is well detected between epochs 
separated by only one month with the Very Long Baseline Array (VLB A) [|l^]. New observations 



increased the time span to 8 years []18[]. The position of Sgr A* relative to the background quasar 
J1745-283 is shown in Figure ^ 

The motion of Sgr A* is 6.379 ± 0.026 mas yr~' and almost entirely in the plane of the 
Galaxy. This gives, converted to Galactic coordinates a proper motion of -6.379 ± 0.026 mas yr~^ 



in Galactic longitude and -0.202 it 0.019 mas yr in Galactic latitude. If one assumes a distance 



to the Galactic Center (Rq) of 7.62 it 0.32 kpc []19|], these motions translate to a velocity of -230 



±10 km s^^ along the Galactic plane and -7.6 ± 0.6 km s^^ out of the plane of the Galaxy. This 
motion can be entirely explained by a combination of a circular rotation of the Local Standard of 
Rest (LSR) ©0 and the deviation of the motion of the Sun from the motion of the LSR. Removing 



the Solar motion relative to the LSR as measured by Hipparcos data in |2C] (5.25 ± 0.62 km s~^ 
in longitude and -7.17 ± km s^^ in latitude) yields an estimate of 0o=225±lO km s^'. Here the 
error is dominated by the uncertainty in the distance Rq to the Galactic Center. The motion of Sgr 
A* out of the plane of the Galaxy is only -0.4 ± 0.8 km s^^ This lower limit can be used to put 
tight constraints on the mass of Sgr A* (for details see [p^]). 

With the accuracy obtainable with VLBI one can in principle measure very accurate proper 
motions for most Local Group members within less than a decade. The main problem so far is 
finding appropriate radio sources. Useful sources would be either compact radio cores or strong 
maser lines associated with star forming regions. Fortunately, in a few galaxies bright masers are 
already known. 

The most suitable candidates for such a VLBI phase-referencing experiment are the strong 
H2O masers in IC 10 (~ 10 Jy peak flux in 0.5 km s"' line) and IC 133 in M33 (~ 2 Jy, the first 
extragalactic maser discovered) [pTl]. The two galaxies belong to the brightest members of the 
Local Group and are thought to be associated with M31. In both cases a relatively bright phase- 
referencing source is known to exist within a degree. In addition, their galactic rotation is well 
known from H I observations. Consequently, M33 and IC 10 seem to be the best known targets for 
attempting to measure Local Group proper motions with the VLBA. 
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3. VLBI observations of M33 and IC 10 

3.1 Observations and Data Reduction 

We observed two regions of H2O maser activity in M33 (M33/19 and IC 133) eight times with 



the NRAO Very Long Basehne Array (VLBA) between March 2001 and June 2005 M33/19 
is located in the south-eastern part of M33, while IC 133 is located in the north-east of M33 (see 
Fig. H). We observed the usually brightest maser in IC 10-SE with the VLBA thirteen times between 



February 2001 and June 2005 |g8|] 



The observations involved rapid switching between the phase calibrator and the target sources 
With source changes every 30 seconds, an integration time of 22 seconds was achieved per scan. 
From the second epoch on, we included geodetic-like observations where we observed for 45 min- 
utes 10-15 strong radio sources (> 200 mJy) with accurate positions (< 1 mas) at different ele- 



vations to estimate an atmospheric zenith delay error in the VLBA calibrator model (see Ill8|, |23|] 
for a discussion). In the second and third epoch we used two blocks of these geodetic observations 
before and after the phase-referencing observations. From the fourth epoch on, we included a third 
geodetic block in the middle of the observation. The data were edited and calibrated using stan- 
dard techniques in the Astronomical Image Processing System (AIPS), and zenith delay corrections 
were applied to improve the accuracy of the phase-referencing. The masers in IC 133, M33/19, and 
IC 10 were imaged with standard techniques in AIPS. 
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Figure 3: Position of two regions of maser activity in M33. Also shown are predicted motions due to rotation 
of the H I disk (shown in blue). The optical image of M33 was provided by Travis Rector (NRAO/AUI/NSF 
and NOAO/AURA/NSF), and the HI image by David Thilker (NRAO/AUI/NSF) and Robert Braun (AS- 
TRON). Taken from 



3.2 Proper Motions of M33/19 and IC 133 

The maser emission in M33/19 and IC 133 is variable on timescales of less than one year. Be- 
tween the epochs, new maser features appeared while others disappeared. However, the motions of 
4 components in M33/19 and 6 components in IC 133 could be followed over all four epochs. The 
component identification was based on the positions and radial velocities of the maser emission. 
Each component was usually detected in several frequency channels. A rectilinear motion was fit 
to each maser feature in each velocity channel separately. Fits with a reduced larger than 3 were 
discarded as they are likely affected by blending. Then, the variance weighted average of all mo- 
tions was calculated. This yields an average motion of the maser components in M33/19 of 35.5 ± 
2.7 /ias yr^' in right ascension and —12.5 ± 6.3 /xas yr^' in declination relative to the background 
source J0137-I-312. For IC 133 one gets an average motion of 4.7 it 3.2 jUas yr^^ in right ascension 
and —14.1 it 6.4 /xas yr^^ in declination. 
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Figure 4: Average positions of the maser in M33/19 in East- West (upper) and North-South (lower) relative 
to the reference source. Taken from [^^. 
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Figure 5: Average positions of the maser in IC 133 in East-West (upper) and North-South (lower) relative 
to the reference source. Taken from 



3.3 Geometric Distance of M33 

The relative motions between M33/19 and IC 133 are independent of the proper motion of 
M33 and any contribution from the motion of the Sun. Since the rotation curve and inclination of 
the galaxy disk are known one can predict the expected relative angular motion of the two masing 



regions. The rotation of the HI gas in M33 has been measured [24] and one can calculate the 
expected transverse velocities of M33/19 and IC 133. This gives a relative motion of 106.4 km s~^ 
in right ascension and 35 km s^' in declination between the two regions of maser activity. 

The radial velocities of the H2O masers in M33/19 and IC 133 and the proximate HI gas are 
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in very good agreement (< 10 km s^^). This strongly suggests that the maser sources are rotating 
with the H I gas in the galaxy. However, while the agreement between the rotation model presented 
in [^] and the radial velocity of the H I gas at the position of IC 133 is also very good (< 5 km s~^), 
there is a difference of ~ 15 km s^^ at the position of M33/19. Hence, we conservatively assume a 
systematic error of 20 km s^ ^ in each velocity component for the relative velocity of the two maser 
components. Comparing the measured angular motion of 30.8 ± 4 /xas yr^' in right ascension with 
the expected linear motion of 106 ± 20 km s^\ one gets a geometric distance of 

D = 730±135±100kpc, 

where the first error indicates the systematic error from the rotation model while the second error 
is the statistical error from the VLBI proper motion measurements. 

After less than three years of observations, the uncertainty in the distance estimate is already 
dominated by the uncertainty of the rotation model of M33. However, this can be improved in 
the near future by determining a better rotation model using higher resolution (e.g.. Very Large 
Array or Westerbork Synthesis Radio Telescope) data of H I gas in the inner parts of the disk. Also, 
the precision of the proper motion measurements will increase with time as f^/^ for evenly spaced 
observations. 

Within the current errors the geometric distance of 730 ib 100 ± 135 kpc is in good agreement 
with recent Cepheid and Tip of the Red Giant Branch (TRGB) distances of 802ib51 kpc [S| and 



794ib23 kpc [26] respectively. 



3.4 Proper Motion of M33 

The observed proper motion Vprop of a maser (M33/19 or IC 133) in M33 can be decomposed 
into three components Vpmp = Vrot + v© + vm33 ■ Here Vrot is the motion of the masers due to the 
internal galactic rotation in M33 and Vq is the apparent motion of M33 caused by the rotation of 
the Sun around the Galactic Center. The last contribution vm33 is the true proper motion of M33 
relative to the Galaxy. 



Since the motion of the Sun Ill8|, gO|] and the rotation of M33 [ ]24[ ] is known, one can calculate 
the two contributions Vrot and v©. Combining these velocity vectors, one gets the true proper motion 
of M33: 



O^M33 — OCprop — OCmt ~ OCq 

= -29.3±7.6^ = -101±35^ 

yr i 

and 

5m33 = 5 prop — 5rot — 5q 

= 45.2±9.1^ = 156±47^. 

yr s 

Finally, the systemic radial velocity of M33 is —179 km s^' [p4|]. The radial component of 
the rotation of the Milky Way toward M33 is —140 ± 9 km s^K Hence, M33 is moving with —39 
± 9 km s^^ toward the Milky Way. This gives now the three dimensional velocity vector of M33 
which is plotted in Fig. ^ The total velocity of M33 relative to the Milky Way is 190 ± 59 km s^^ 
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Figure 6: The position of the maser in IC 10 relative to a continuum background source in East- West (red 
triangles) and North-South (blue circles). The lines show a variance weighted linear fit to the data points. 
Taken from 



3.5 Proper Motion of IC 10 



In IC 10 only the strongest maser component was detected in all epochs. The position offsets 



of this maser feature in IC 10 are shown in Fig. |3.5| . The uncertainties in the observations of the first 
epoch are larger than the others, because no geodetic-like observations were done to compensate 
the zenith delay errors. A rectilinear motion was fit to the data and yielded a value of 6±5 /xas yr^' 
toward the East and 23 ±5 jxas, yr^^ toward the North. 

Once again, the contributions v,.„; and Vq can be calculated from the known motion of the Sun 



[ ]1 8| , gOp and the known rotation of IC 10 [ |29| , pO| ]. The true proper motion of IC 10 is then given 
by: 



OI/CIO — OCprop — OCrot — OCq 

= -39±9Atasyr"' = -122 ±31 km s^^ 

and 

5/ClO = Spi-op — Srot — Sq 

= 31±8^as yr"^ =97±27kms~^ 



The measured systematic heliocentric velocity of IC 10 (—344 it 3 km s \ [ 3 1 ]) is the sum of 
the radial motion of IC 10 toward the Sun and the component of the solar motion about the Galactic 
Center toward IC 10 which is -196±10 km s-\ Hence, IC 10 is moving with 148 it 10 km s 
toward the Sun. The proper motion and the radial velocity combined give the three-dimensional 
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Figure 7: Schematic view of the Local Group with the space velocity of M33 and the radial velocity of 
Andromeda. The blue cross marks the position of the Local Group Bary center (LG BC) according to [p7[]. 
Taken from [||]. 
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space velocity of IC 10. This velocity vector is shown in the schematic view of the Local Group in 
Fig. The total velocity is 215±42 km s^^ relative to the Milky Way. 



4. Local Group Dynamics and Mass of M31 

If IC 10 or M33 are bound to M3I, then the velocity of the two galaxies relative to M31 must 
be smaller than the escape velocity and one can deduce a lower limit on the mass of M3 1 : 



Mm31 > 



2G 

A relative velocity of 147 km s^' - for a zero tangential motion of M31 - and a distance of 262 
kpc between IC 10 and M31 gives a lower hmit of 6.6 x IO'^Mq. One can repeat this calculation 
for any tangential motion of M31. The results are shown in Fig. ^ (top). The lowest value of 0.7 
X IO^'Mq is found for a tangential motion of M31 of -130 km s^' toward the East and 35 km s^' 
toward the North. 

For a relative motion of 230 km s^' between M33 and M31 - again for a zero tangential 



motion of M31 - and a distance of 202 kpc, one gets a lower limit of 1.2 x lO'^M© |22]. Fig. || 
(top) shows also the lower hmit of the mass of M3 1 for different tangential motions of M3 1 if M33 
is bound to M31. The lowest value is 4 xlO^^M© for a tangential motion of M31 of -115 km s^^ 
toward the East and 160 km s^^ toward the North. 



In [ |32| ] it is found that proper motions of M31 in negative right ascension and positive dec- 
lination would have lead to close interactions between M3 1 and M33 in the past. These proper 
motions of M31 can be ruled out, since the stellar disk of M33 does not show any signs of strong 
interactions. 

Thus, we can rule out certain regions in Fig. |[ This results in a lower limit of 7.5 x lO^^M© for 
M31 and agrees with a recent estimate of 12. 3;^^^^ x 10'^ M0 derived from the three-dimensional 



positions and radial velocities of its satellite galaxies [33]. 



5. Summary 

More than 80 years after van Maanen's observation, we have succeeded in measuring the 
rotation and proper motion of M33 as well as the proper motion of IC 10. These measurements 
provide a new handle on dynamical models for the Local Group and the mass and dark matter halo 
of Andromeda and the Milky Way. 

We have presented astrometric VLBA observations of the H2O masers in the Local Group 
galaxies M33 and IC 10. We measured the proper motion of the masers relative to background 
quasars. Correcting for the internal rotation of M33, IC 10 and the rotation of the Milky Way these 
measurements yield proper motions of the two galaxies. The total space velocities relative to the 
Milky Way of M33 and IC 10 are 190±59 km s^^ and 215±42 km s^^, respectively. If IC 10 and 
M33 are bound to M31, one can calculate a lower limit of the mass of M31 of 7.5 x 10^ ^M0. 

Further VLBI observations within the next few years and an improved rotation model have the 
potential to improve the accuracy of the distance estimate to less than 10%. At least one additional 
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Figure 8: Top: Lower limit on the mass of M3 1 for different tangential motions of M3 1 assuming that M33 
(dashed) or IC 10 (soHd) are bound to M31. The lower Umits are (4, 5, 7.5, 10, 15, 25)x lO'^M© for M33, 
and (0.7, 1, 2.5, 5, 7.5, 10, 15, 25)x1O"M0 for IC 10, rising from inside. The colour scale indicates the 
maximum of both values. Bottom: The colour scale is the same as above and gives the lower limit on the 
mass of M3 1 . The contours show ranges of proper motions that would have lead to a large amount of stars 
stripped from the disk of M33 through interactions with M31 or the Milky Way in the past. The contours 
delineate 20% and 50% of the total number of stars stripped [Q. These regions can be excluded, since the 
stellar disk of M33 shows no signs of such interactions. Taken from 
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maser source exists in M33 [21] that will be used in the future to increase the accuracy of the 
measurements. A third region of maser activity will also help to check for non-circular velocities 
of the masers. 

Unfortunately, no maser emission could be found in Andromeda despite intensive searches. 
Hence, the proper motion vector of Andromeda is still unknown. In the near future, new technical 
developments using higher bandwidths will increase the sensitivity of VLBl. This will allow the 
detection of radio emission from the central black hole of Andromeda, M31*, with flux densities 
of ~30 /iJy [p4|], to measure the proper motion of Andromeda. 

Today we are only able to study the extreme (bright) tip of the maser luminosity distribution 
for interstellar masers. The Square Kilometer Array (SKA), with substantial collecting area on 
intercontinental baselines and a frequency coverage up to 22 GHz [^], will provide the necessary 
sensitivity to detect and measure the proper motions of a much greater number of masers in active 
star forming regions in the Local Group. 
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